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ABSTRACT 



We report on the detection of a 400.99018734(1) Hz coherent signal in the Rossi X-ray Timing Explorer light curves of the recently 
discovered X-ray transient, IGR J17498-2921. By analysing the frequency modulation caused by the orbital motion observed between 
August 13 and September 8, 201 1, we derive an orbital solution for the binary system with a period of 3.8432275(3) hr. The measured 
mass function, f(M 2 , Mi, i) = 0.00203807(8) M G , allows to set a lower limit of 0.17 M G on the mass of the companion star, while an 
upper limit of 0.48 M Q is set by imposing that the companion star does not overfill its Roche lobe. We observe a marginally significant 
evolution of the signal frequency at an average rate of -(6.3 ± 1.9) x 10~ 14 Hz s -1 . The low statistical significance of this measurement 
and the possible presence of timing noise hampers a firm detection of any evolution of the neutron star spin. We also present an 
analysis of the spectral properties of IGR J17498-2921 based on the observations performed by the Swift-X-ray Telescope and the 
i?AT£-Proportional Counter Array between August 12 and September 22, 201 1. During most of the outburst, the spectra are modeled 
by a power-law with an index T w 1.7-2, while values of » 3 are observed as the source fades into quiescence. 

Key words, stars: neutron — stars: rotation — pulsars: individual (IGR J17498-2921) — X-rays: binaries 



1. Introduction 

The discovery of a 401 Hz coherent signal from the X-ray tran- 
sient SAX J1808.4-3658 dWiinands & van der KlislfT998h con- 
firmed that old neutron stars (NS) in low mass X-ray binaries 
(LMXB) are rapidly rotating objects, spun-up by the accretion 
of angular momentum from matter in-falling through a disc. 
Sources showing a coherent signal with a period of few ms are 
called accreting millisecond pulsars (AMSP); they are all rela- 
tively faint X-ray transients attaining peak X-ray luminosities of 
a few x 10 36 erg s , with outbursts lasting up to a few months. 

Here we report on the discovery of 401 Hz pulsations from 
the X-ray transient IGR J 17498-2921, which was first de- 
tecte d on 2011 Aug 11 9 thanks to an INTEGRAL observa- 
tion (Gib aud et al.ll201 U the epochs in this paper are referred 
to the barycentric dynamical time, TDB). A subsequent Rossi 
X-ray Timing Explorer (RXTE) observati on allowed the dis - 
covery of a coherent signal at 401 Hz (Papitto et all 1201 lbl) . 
bringing to 12 the number of known AMSPsQ. The source po- 
sition was det ermined by follow-up S wift (Bo zzo et al.l 1201 ll) 
and Chandra (Chakrab artv et al.l 1201 lb observations. The lat- 
ter estimate, RA=17 /! 49 m 55 s .35, DEC=-29° 19' 19.6", with 
an uncertainty of 0.6" at a 90% confide nce level, is consid- 
ered here. Bursts were detected bv lFerrigno et al.l (120 11) during 
INTEGRAL observations, while iLinares et alj d201 lal) found 
oscillations at a frequency consistent with the spin frequency 



1 T he properties of the inter mittent pulsars, Aql X- l iCasella et alJ 
2008) and SAX J1748.9-2021 jAltamirano et alj|2008l) , seem distinct 
enough not to count them in this number. 



during a burst observed by RXTE, as well as a photospheric ra- 
dius expansion episode, from which a distance estimate of 7.6 
kpc was derived. The transient retur ned to quiescence on 2011 
Sep 19, after a 37 day-long outburst dLinares et al.ll2.Ql lbl) . 

The analysis of the RXTE observ ations performed during 
the first few days of the outburst allowed Markwardt e t all (120 1 ll) 
to prop ose a set of candidate orbital solutions, and Pamtto et al. 
(1201 lal) to give the first preliminary solution of a 3.84 hr orbit 
of the NS. Here we present the first analysis of the IGR J 17498- 
2921 observations performed by RXTE and Swift, focusing on 
the properties of the 401 Hz signal observed by RXTE and de- 
riving a full orbital and timing solution for the pulsar. 



2. Observations and data analysis 

We processe d data obtained by th e Swift/X-ray Telescope (XRT; 
0.2-10 keV. lGehrels etal.112004 between 2011 Aug 12.6 and 
Sep 22.7 (MJD 55785.6-5 5826.7), for an expos ure of 48.3 ks, 
using standard procedures dBurrows et al.l 12005) and the latest 
calibration files available at the time of the analysis (September 
2011). We considered observations performed in photon count- 
ing mode, retaining the full spectroscopic and imaging capabil- 
ities of the XRT with a time resolution of 2.5 s, and windowed 
timing mode, which preserves only one imaging axis to allow a 
read out every 1.7 ms. According to the recommendations of the 
XRT calibration tearrQ, we produced 1-10 keV energy spectra 
grouping channels to contain at least 50 counts each, adding an 
estimated systematic error of 3% to each bin. 



see http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/swift/docs/xrt/SWIFT-X 
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Reduction and analysis of the data obtained by the 
RXTE/Proportiona l Counter Array (PCA; 2-60 keV; 
Uahoda et all 120061) between 2011 Aug 13.1 to Sep 22.4 
(MJD 55786.1-55826.4), for an exposure of 94.2 ks, was 
carried out with the tools available in HEASOFT (ver. 6.11). 
Spectra were extracted from data encoded in Standard 2 mode 
(time resolution of 16 s and 128 spectral channels). We retained 
only photons collected by the top layer of the Proportional 
Counter Unit (PCU) 2 in the 2.5-30 keV band, adding a 
systematic error of 0.5% to each spectral channel, to get the 
most reliable description of the PCA energy redistribution^! The 
background was estimated with the 'bright' model. A timing 
analysis was performed on 2-60 keV photons observed by 
all the PCUs switched on during each of the observations, to 
achieve the highest possible quality counting statistics. Data 
processed in fast encoding modes such as event, single bit 
(both with 122 /is time resolution), and good xenon (1 /js time 
resolution) are considered. 

Since in these paper we deal with the source properties dur- 
ing its persistent (i.e., non-bursting) emission, we discard 5s of 
data before each of the seven bursts detected by RXTE, as well 
as before the two observed by Swift, and an interval of variable 
length depending on the burst properties (~ 60 s), after. 

2. 1 . Light curves and spectra 

The 1-10 keV Swift/XRT spectra are satisfactorily modelled 
by an absorbed power law. The best-fit value of the absorption 
column ranges from 2.5 to 3.5xl0 22 crrT 2 . Since these values 
are compatible and much less acc urate than the Chandra esti- 
mate given by iTorres et al.l (1201 lb , nH = (2.87 + 0.04) x 10 22 
cm -2 , we fixed the absorption column to this value in the fitting. 
The 2-10 keV absorbed flux is plotted as red circles in panel 
(a) of Fig. [TJ together with a dashed line representing the av- 
erage flux trend obtained by fitting with a sixth-order polyno- 
mial also the flux observed by the RXTE/PCA in the same en- 
ergy band (blue triangles; see below). We measured an excess 
of 2.1(2) x 10~ 10 erg cirT 2 s with respect to this trend dur- 
ing the Swift observation starting on 2011 Aug 20.5, a value 
significantly larger than the average RMS scatter in the flux ob- 
served, 3 .5 x 1 0~ 1 1 erg cirT 2 s~ 1 . The flux during this observation 
of F™ A = (5.8 + 0.2) x 10" 10 erg cirT 2 s -1 , is the peak value 
observed during the outburst and corresponds to an unabsorbed 



(5.1 + 0.2) x 10 36 d 2 erg s -1 , where d 



7.6 



luminosity of ^-ui 
is the distance to the source in units of 7.6 kpc. An excess of 
1.0(2) x 10~ 10 erg cirT 2 s is also observed during the obser- 
vation starting on 2011 Aug 21.5. By comparing these observa- 
tions with the nearest ones, a timescale as low as « 1000 s is 
found for the flux variability. Spectral variability is not observed 
as the photon index, F, has values that are consistent with those 
observed during most of the Swift/XRT observations (between 
1.7 and 2; see red circles in panel (b) of Fig. [TJ). A steepening 
of the power law is instead observed as the transient fades to 
quiescence. The source was no longer detected after Sep 19; by 
co-adding the observations performed from Sep 19.4 to 22.7, we 
obtained an upper limit of 3.2 x 10~ 13 erg cirT 2 s to the 2-10 
keV flux (3 cr c. 1., assuming a power-law with index F =1.8). 

The X-ray transient IGR J17498-2921 lies very close to the 
central regions of the Galactic bulge (1=0.1559°, b=-1.0038°) 
and the emission of the Galactic ridge and numerous nearby 
sources contaminate the observations in the field of view of the 
PCA, which is about 1 by 1 degrees with a hexagonal shape. 
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Fig. 1. Panel (a): absorbed 2-10 keV flux of IGR J17498-2921 
as observed by the Swift/XRT (red circles) and by RXTE/PCA 
(blue triangles), and its best-fit sixth-order polynomial relation 
(dashed line). The PCA fluxes were obtained after subtracting 
the spectrum observed during the interval 2011 Sep 20.5-22.4 
(days 39.5 and 41.4 in the scale used here), when the source 
already faded into quiescence; panel (b): photon index of the 
power-law spectra observed by Swift/XRT (red circles) and 
RXTE/PCA (blue triangles); panel (c): RMS amplitude of the 
401 Hz signal as observed by the PCA; panel (d): pulse phases 
obtained by folding around v = 400.99018734 Hz the PCA 
time-series corrected for the orbital motion by using the param- 
eters listed in Table [TJ together with the best-fit quadratic model 
plotted as a dashed line. 



iRevnivtsev et al.l(| 2009) estimated a 2-10 keV surface brightness 
of (8.6 + 0.5) x 10" 11 erg cirT 2 s _1 deg~ 2 from a nearby field not 
containing bright sources. Moreover, ten sources with 2-12 keV 
flux exceeding 5 x 10" 11 erg cirT 2 s _I were found in a region of 
1° around the position of IGR J 17498- 2921 in the fifth X MM- 
Newton serendipitous source catalogue dWatson et al.1 2009); the 
brightest of these sources is the persistent burster 1A 1742-294 
with a 2-10 keV flux between 2 and 4xl0~ 10 erg cirT 2 s _1 . To 
estimate the contamination of the PCA spectra, we considered 
observations performed between 2011 Sep 20.5 and Sep 22.4, 
for an exposure of 3.7 ks, when the transient had already faded 
into quiescence (see above). The net count-rate observed by the 
top layer of the PCU2 during these observations was 32.5(2) s _1 
(36.5(2) s _1 , if all the layers are considered) and the absorbed 



see http://www.universe.nasa.gov/xrays/programs/rxte/pca/doc/rmf/ 



2-10 keV flux is (3.5 + 0.1) x 10~ 10 erg cirT 2 s _1 . We consid- 
ered the spectrum obtained during these PCA observations as an 
additional background to those performed previously. After this 
subtraction, all the PCA spectra could be accurately described 
by an absorbed power-law, with the value of nH fixed to the 
Chandra estimate. The 2-10 keV fluxes and the best-fit values of 
the power-law photon index, ranging from 1.8 to 2, are plotted 
as blue triangles in panel (a) and (b) of Fig.Q] respectively. We 
found that the best-fit models of the simultaneous Swift -XRT 
and RXTE-PCA spectra are consistent with each other. 

The bolometric fluence of the outburst was estimated to be 
T = (4.1 ± 0.6) X 10~ 3 d 2 6 erg cm -2 by integrating the aver- 
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age trend of the 2-10 keV source flux (dashed line in panel (a) 
of Fig.Q} over the whole outburst length, and scaling the value 
obtained by a bolometric correction factor (cb„i = 6.0) equal 
to the average observed ratio of the unabsorbed flux extrapo- 
lated to the 0.05-100 keV band to the absorbed flux observed 
in the 2-10 keV band. Assuming that AM = R„T/GM t , the 
mass accreted by the NS during the outburst is estimated to be 
AM = (8±l)xl0~ n d 2 6 M for aM* = 1.4M Q NS wi th a radius 
of R* = 10 km, a value typical of outbursts of AMSPs ( Gal lowavl 
120061) . However, these values can be overestimated by a factor of 
~ 2 because of the unknown spectral distribution of the source 
outside the energy band considered here. Describing the broad- 
band spectrum with a Comptonization model with input photons 
and hot scattering electrons with temperatures of 0.5 and 50 keV, 
respectively, would rather yield a value of Cf,oi — 3. 

2.2. Timing analysis 

To perform a timing analysis of the 401 Hz signal, we reported 
the photon arrival times to the barycentre of the Solar System 
by using the Chandra position dChakrabartv et alJl201 lh and ap- 
plied fine clock corrections using the latest calibration file issued 
by the PC A team at the time of the analysis (September 201 lfl 
to reach an accuracy in absolute timing of between 5 and 8 /is. 

A signal at a frequency of between 400.95 and 401.0 Hz is 
observed in power spectra produced over 256 s intervals, at a 
Lehay-normalised power of « 60-70. A first solution of the or- 
bital motion of the NS was obtained by performing an epoch- 
folding search over 256 s intervals and modelling the variations 
in the signal period caused by the orbital motion. We thus ob- 
tained a first estimate of the orbital parameters of the source, 
such as the orbital period, P ,b =13835.65(5) s, the projected 
semi-major axis of the NS orbit, asim'/c =0.3651(2) lt-s, the 
epoch of zero mean anomaly, T*, and upper limits to e sin to and 
e cos u>, where e is the eccentricity and u is the longitude of the 
periastron measured from the ascending node, as well as of the 
barycentric spin frequency, vo = 400.99020(2) Hz. 

The times of arrival of X-ray photons were then corrected by 
using this first orbital solution, and the time series folded into 16 
phase bins around the best estimate of the spin period, over 1500 
s time intervals. The presence and significance o f pulsations 
in the folded light curves were assessed following Leah y et al.l 
(119831) . The background-subtracted 2-60 keV RMS amplitude 
of the observed profiles is plotted in panel (c) of Fig. Q] To take 
into account the contribution of the the Galactic ridge and of 
nearby sources to the PCA count rate, we considered an addi- 
tional background level of 36.5(1) s (determined from PCU2 
observations of the field, see Sec. 12. Il l for every active PCU. This 
is only an approximation since the responses and background 
levels of the various PCUs differ by about 10 per cent, intro- 
ducing an error in the estimated amplitude that is comparable to 
the statistical uncertainties in the various measurements. The in- 
trinsic pulse RMS amplitude thus evaluated takes values ranging 
from 0.06 to 0.1 1, with relative uncertainties of the order of 10- 
15 per cent. The pulse profiles can generally be modelled by a 
sinusoid. A second harmonic with an RMS amplitude between 1 
and 2 per cent, and a relative uncertainty of * 25 per cent, is de- 
tected only during the brightest observations. Pulsations are not 
detected during the observations performed on 201 1 Sep 14 and 
15 (MJD 558 18-558 19), with an upper limit of about 20 per cent 
to the pulsed fraction, reflecting our approximate knowledge of 
the true background level, especially when the source flux is low. 

4 see http://heasarc.gsfc.nasa.gov/docs/xte/time_news.html 



Table 1. Spin and orbital parameters of IGR J17498-2921. 



v (Hz) 


A A A finn 1 OTi A i Iwl A— H i O * / 1 A— X 

400.99018/34 ± 1 X 1U ± o X 11) 


<v> (Hz s -1 ) 


(-6.3 ± 1.1 + 1.5) X 10~ 14 


a sin i/c (lt-s) 


0.365165(5) 


Porb (S) 


13835.619(1) 


r* (MJD) 


55785.0600534(8) 


e 


< 8 x 10" 5 


f(m 2 ,m u i) (M B ) 


0.00203807(8) 


^/d.o.f. 


76/50 



Notes. The reference epoch for the solution is f = MJD 55786.124. 
Numbers in parentheses are l<x errors in the last significant digit, while 
upper limits are evaluated at the 3<r confidence level. The uncertainties 
reported for vo and < v > are the statistical and the systematic uncer- 
tainty caused by the position error, respectively. 

The estimates of the spin and orbital parameters of the pulsar 
were r efined by applying timing techniques to the pulse phase 
delays dDeeter et al.Hl 98 ll) . i.e., modelling their evolution as 

</>(t) - 0(f) = [vo - v F ](t - T ) + - <v> (f - T ) 2 + R orb {t), (1) 

where T = MJD 55786.124 is the reference epoch of the fold- 
ing, Vf is the folding frequency, < v > is the average curva- 
ture of the trend followed by the phases, equal to the average 
spin frequency derivative under the assumption that the spin is 
well-tracked by the pulse phases, and R or b(t) describes a phase 
modulation caused by the slight differences between the actual 
orbital parameters of the system and those used to correct the 
time series, i.e., P or b, a sin i/c, T*, e sin a>, and e cos to. The best- 
fit parameters that we found are listed in Table [1] while the 
pulse phases obtained by folding the RXTE light curves around 
v = 400.99018734 Hz are plotted in panel (d) of Fig. [Q to- 
gether with the best-fit model. The modeling of the pulse phase 
delays of IGR J 17498-2921 with Eq. (\} is satisfactory since 
;f 2 /d.o.f. = 76/50 - 1.5. We measured a statistically significant 
curvature of the phase evolution < v >= (-6.3 + 1.1) x 10~ 14 
Hz s _1 , which is also indicated by the improvement in the chi- 
squared of the model obtained by adding the quadratic term to 
a linear model, Ay 2 = -47. The interpretation of this curva- 
ture in terms of a spin frequency derivative should indeed be 
taken with caution. If the phases relative to the first RXTE ob- 
servation (spanning 201 1 Aug 14.124-14.145; MJD 55786.124- 
55786.145) were removed from the fit, the estimate of < v > 
would become much less significant, (-3.9 + 1.0) x 10~ 14 Hz 
s . Moreover, the uncertainty in the source position introduces 
a systematic error of the same order to this es timate. By consid- 
ering the uncertainty in the position given by IChakrabartv et al.l 
(1201 ll) thanks to an observation perfomed by the HRC-S detec- 
tor on-board Chandra (0.6" at 90% confidence level, which cor- 
responds to 0.4" at 1 cr c.l. according to the distribution of aspect 
offset of the detectotQ), we derived a systematic error in the esti- 
mates of the spin frequency and of its der ivative equal to 8 X 10~ 8 
Hz an d 1.5 x 10~ 14 Hz s _1 , respectively dLvne & Graham-S mith 
1990). Adding in quadrature this error to the statistical uncer- 
tainty in < v>, we found that the detection of a frequency deriva- 
tive is only marginally significant. The modulus of the value ob- 
tained by fitting the whole data-set, < v>= (-6.3 ± 1.9) x 10~ 14 
Hz s , is 3.3 times larger than its overall uncertainty, while the 
estimate obtained by removing the first RXTE observation from 
the sample becomes compatible with zero within the error. 

5 see http://cxc.harvard.edu/cal/ASPECT/celmon/ 
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3. Discussion and conclusions 

We have presented an analysis of the observations performed by 
RXTE and Swift during the outburst shown between August and 
September 2011 by the newly discovered AMSP IGR J17498- 
2921. We have reported the discovery of a coherent signal at a 
frequency of 401 Hz in the RXTE light curves, as well as derived 
an accurate orbital solution for the 3.84 hr binary system. 

The pulse profile shown by IGR J 1 7498-292 1 is modelled by 
a sinusoid with an RMS amplitude between 6 and 1 1 per cent. 
A second harmonic with an amplitude about ten times lower 
is detected only in a subset of observations. Nearly sinusoidal 
profiles, with a similar ratio of the second to the first harmonic 
amplitudes, had alread y been observed from XTE J175 1—305 
(Mar kwardt etaDl2002h an d IGR J0029 1+5934 dGallowav et alJ 
2005; Falanga et al. l2005h . The pulse phases shown by IGR 
J 17498-2921 are quite stable and their evolution can be closely 
described by a polynomial of low (< 2) order, in addition 
to the delays introduced by the orbital motion. Interestingly, 
this property is shared b y IGR J 17498-2921 wit h XTE J 1751- 
305 dPapitto et al.ll2008l) and IGR J00291 +5934 (iFalanga et al l 
20051: iBurderi et al 1 120071 lPatrunoll20lol: lHartman et all 1201 U 
Papit to et alj|201 lq) . while the phases of other AMSPs are in- 
stead affected by timing noise (see, e.g.. lHartman et alJ [2008: 
Rigg io etalJ 120081) . We note how the phases of IGR J17498- 
2921 are affected by an error of «0.01, a value consistent with 
the phase uncertainty of a noisy pulsar suc h as SAX J 1808.4- 
3658 (see Table 2 of lHartman et all 120081) . such that the en- 
hanced phase stability observed in this case does not seem to 
be caused (at least entirely) by a lower quality pulse statistics. 
The evolution of the pulse phases of IGR J 17498-2921 is best- 
fitted by a parabola with positive curvature, indicating a NS spin- 
down at a rate of (-6.3 ± 1.9) x 10 -14 Hz s _I if interpreted in 
terms of the NS spin frequency derivativ e. A similar trend was 
already observed f rom other two AMSPs (jGallowav et al. 2002; 
Papit to et alJ2007h . but the low significance of this measurement 
and the possible presence of residual timing noise prevents us 
from drawing a firm detection of a NS spin evolution. 

The value of the pulsar mass function /(M2, Mi , ;') 2x 10~ 3 
M allows us to set constraints on the nature of the compan- 
ion star. Since no eclipses are observed during RXTEfPCA and 
Swift/XRT observations performed when the NS is behind the 
companion (i.e., at true orbital longitudes ~ 90°), we have in- 
ferred that the inclination of the system is less than 80°. This 
provides a lower limit to the companion mass of mi i 0.066 + 
0.077m 1, where m\ and m 2 are the masses of the two stars in 
solar masses (m 2 i 0.17 for m\ = 1.4). An upper limit to the 
companion mass is obtained by assuming that it does not over- 
fill its Ro che lobe, t he siz e of which is estimated by the re- 
lation of iPaczvnskH dl97ll) . R L2 = QA62a[m 2 l{m x + m 2 )] 1/3 , 
where a is the binary separation. By using the third Kepler law 
to relate the size of the orbit to both the total mass of the sys- 
tem and the observed orbital period, one obtains R2 i Rli = 
0.573 R Q m)P ■ This relation is plotted in Fig.|2]together with the 
zero age main sequence (ZAMS) mass-radius relation given by 
IChabrier & Baraffd (l2000h . By assuming that the companion star 
fills its Roche Lobe and follows the ZAMS mass-radius relation, 
we derived values of m 2 = 0.48 and i = 24.6° for the companion 
star mass and the system inclination, respectively. While heavier 
companion stars can be excluded, lower values are possible if 
the companion is bloated with respect to its thermal equilibrium 
radius. If a NS with a mass of 2 M Q is considered, the lower 
limit to the system inclination set under these assumptions in- 
creases to 28.8°. Among the scenarios explaining the evolution 
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Fig. 2. Size of the Roche lobe of the donor of IGR J 17498- 
2921 (solid curve) and ZAMS mass-radius relation given by 
IChabrier & Ba raffe (2000, dotted line). The masses and radii of 
donor stars filling their Roche lobe, calculated from the mea- 
sured mass function assuming m\ = 1.4 and a set of values of 
inclination, are indicated as black points. The left shaded region 
is excluded by the observed mass function and by the absence of 
eclipses (see text), while the right one is excluded by the condi- 
tion that the companion star does not overfill its Roche lobe. 



of LMXBs (see Deloye 2008, for a review), the orbital period of 
IGR J17498-2921 and the limits to the mass of the companion 
derived here are consistent with a cataclysmic variable-like evo- 
lutionary path, with the low-mass donor that made contact with 
its Roche lobe before evolving to t he red giant branch (see, e.g., 
the sequences plotted in Fig. 2 of Podsiadl owski et al.ll2002l) . 
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